Finite Volume Methods and Adaptive Refinement for Global Tsunami Propagation and Local Inundation
نویسندگان
چکیده
The shallowwater equations are a commonly accepted approximation governing tsunami propagation. Numerically capturing certain features of local tsunami inundation requires solving these equations in their physically relevant conservative form, as integral conservation laws for depth and momentum. This form of the equations presents challenges when trying to numerically model global tsunami propagation, so often the best numerical methods for the local inundation regime are not suitable for the global propagation regime. The different regimes of tsunami flow belong to different spatial scales as well, and require correspondingly different grid resolutions. The long wavelength of deep ocean tsunamis requires a large global scale computing domain, yet near the shore the propagating energy is compressed and focused by bathymetry in unpredictable ways. This can lead to large variations in energy and run-up even over small localized regions. We have developed a finite volume method to deal with the diverse flow regimes of tsunamis. These methods are well suited for the inundation regime—they are robust in the presence of bores and steep gradients, or drying regions, and can capture the inundating shoreline and run-up features. Additionally, these methods are well-balanced, meaning that they can appropriately model global propagation. To deal with the disparate spatial scales, we have used adaptive refinement algorithms originally developed for gas dynamics, where often steep variation is highly localized at a given time, but moves throughout the domain. These algorithms allow evolving Cartesian sub-grids that can move with the propagating waves and highly resolve local inundation of impacted areas in a single global scale computation. Because the dry regions are part of the computing domain, simple rectangular cartesian grids eliminate the need for complex shoreline-fitted mesh generation. Science of Tsunami Hazards, Vol. 24, No. 5, page 319 (2006)
منابع مشابه
High Resolution Methods and Adaptive Refinement for Tsunami Propagation and Inundation
We describe the extension of high resolution finite volume methods and adaptive refinement for the shallow water equations in the context of tsunami modeling. Godunov-type methods have been used extensively for modeling the shallow water equations in many contexts, however, tsunami modeling presents some unique challenges that must be overcome. We describe some of the specific difficulties asso...
متن کاملA new numerical model for simulating the propagation of and inundation by tsunami waves A new numerical model for simulating the propagation of and inundation by tsunami waves
A new unstructured grid, finite volume ocean model is described, that is suitable for flooding and drying problems. Here we derive the finite volume analogue of the P 1 −P1 finite element, by interpreting the advective term in the continuity equation in a flux sense. A corresponding non-overlapping control volume is then selected for the momentum equation. The resulting model employs the median...
متن کاملNumerical Modeling of Tsunami Waves Associated With Worst Earthquake Scenarios of the Makran Subduction Zone in the Jask Port, Iran
The recent studies show that the past researches may have significantly underestimated earthquake and tsunami hazard in the Makran Subduction Zone (MSZ) and this region is potentially capable of producing major earthquakes. In this study, the worst case possible earthquake scenarios of the MSZ are simulated using fully nonlinear boussinesq model to investigate tsunami hazards on the Jask Port, ...
متن کاملThe VOLNA code for the numerical modelling of tsunami waves: generation, propagation and inundation
A novel tool for tsunami wave modelling is presented. This tool has the potential of being used for operational purposes: indeed, the numerical code VOLNA is able to handle the complete life-cycle of a tsunami (generation, propagation and run-up along the coast). The algorithm works on unstructured triangular meshes and thus can be run in arbitrary complex domains. This paper contains the detai...
متن کاملThree Dimensional Numerical Simulation of Tsunami Generation and Propagation Due to Makran Subduction and run-up on Chabahar Bay and Makran Coasts
Makran subduction located at the northwest of the Indian Ocean nearby the southern coast of Iran and Pakistan. Makran subduction is the source of tsunamis that threaten southern coast of Iran. In this article, generation and propagation of 1945’s tsunami initiated by Makran subduction is simulated. For the three dimensional generation of the wave, advanced algorithm of Okada is adopted. The CFD...
متن کامل